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Abstract. We show analytically, numerically, and experimentally that a transversely stable 
one- dimensional [(1+1)D] bright Kerr soliton can exist in a 3D bulk medium. The 
transverse instability of the soliton is completely eliminated if it is made sufficiently 
incoherent along the transverse dimension. We derive a criterion for the threshold of 
transverse instability that links the nonlinearity to the largest transverse correlation 
distance for which the ID soliton is stable. 
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Research on optical spatial solitons has made much progress during the past decade: 
many new systems that support solitons have been identified, solitons of more than one 
transverse dimension have been demonstrated, and a whole range of soliton interactions was 
explored [1]. Despite the diversity of the physical systems that support them, solitons are a 
universal phenomenon and share many common features [1], one of which is Transverse 
Instability (TI) [2-10]. TI is a symmetry breaking instability: almost all solitons [11] of a 
particular dimension that propagate in a higher dimension system (by having a uniform 
wavefunction in the additional dimension/s) are unstable to perturbations in the dimension/s 
in which they are uniform. TI occurs because perturbations in the dimension of uniformity 
have nothing to restrain them from growing (driven by the nonlinearity) and breaking the 
soliton up. 

In the particular case of a spatial optical (1+1)D soliton that is self- trapped in one 
dimension, x, is uniform in the transverse dimension, y, and is propagating along z, TI causes 
the soliton to break up along y into an array of 2D filaments [2-10]. The 
transverse-wavelength of these perturbations is usually much larger than the soliton width 
[2-4]. Transverse instability is especially severe for Kerr nonlinearities, and prohibits spatial 
ID Kerr solitons in a bulk medium. This is why spatial (1+1)D Kerr solitons have to be 
launched in a planar waveguide configuration, in which the ^-confinement is much narrower 
than the self- trapped (soliton) width in x [5,6]. TI actually occurs for solitons in any 
nonlinearity, including, for example, quadratic solitons [7] and photorefractive solitons [8,9]. 
Interestingly, saturation arrests transverse instability [10] but never completely eliminates it. 
In fact, it is the suppression of TI due to saturation that facilitates the observation of stable 
ID solitons in a bulk photorefractive crystal for more than 10 diffraction lengths [8]. 

Thus far, in order to avoid TI, experiments with ID solitons were either conducted in 
planar waveguides [5,6] or in nonlinearities in which TI was greatly suppressed [8]. Here we 
demonstrate how to produce a truly stable stripe Kerr soliton propagating in a 3D bulk 
medium without suffering from transverse instability. We show that if the soliton is made 
"sufficiently" incoherent in its transverse dimension, y, then TI is completely eliminated. 
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Before proceeding to TI, it is essential to recall incoherent solitons. These are solitons 
made of partially incoherent light, or self- trapped weakly- correlated wave-packets [12]. They 
are multimode (speckled) beams of which the instantaneous intensity distribution is varying 
randomly with time. If such beams are launched into a non-instantaneous self-focusing 
medium, such that the response time of the nonlinearity greatly exceeds the characteristic 
fluctuation time, then self-focusing is driven solely by the average intensity. In this manner, 
the incoherent beam induces a multimode waveguide and traps (guides) itself in it by properly 
populating the guided modes, thus forming an incoherent soliton [12-18]. 

A clue that TI could be completely eliminated for solitons was given by two recent 
discoveries: Modulation Instability (MI) of incoherent light [19] and elliptical incoherent 
solitons [14,20]. MI belongs to the same family of symmetry breaking instabilities as TI 
does, and it occurs when a plane wave (or a very broad beam or pulse) is launched into a 
self-focusing medium. It is well established that if the plane wave is fully coherent, it breaks 
up into a train of filaments due to MI. Recently, it has been shown theoretically and 
experimentally [20], that MI does exists also for incoherent light, but it occurs only if the 
nonlinearity exceeds a well-defined threshold. The MI threshold is determined by the 
coherence of the light. If the nonlinearity is below threshold, then MI is eliminated and the 
wave is stable. This generic idea has enabled the observation of anti-dark solitons [21], which 
were thought to be unstable in conservative nonlinear systems [22]. The new finding of 
incoherent elliptical solitons, which can exist in logarithmic [14] and saturable [20] 
nonlinearities, is based on having different coherence properties for the two transverse 
dimensions of self-trapping (anisotropic coherence function [14,20]). Combing these ideas, 
one can generate a ID soliton that is fully coherent in x (direction of trapping), partially 
incoherent but uniform in y, and propagating along z. The intimate relation between TI and 
MI suggests that TI of incoherent beams should also exhibit a threshold for its existence. 
Therefore, if the degree of coherence my is such that TI is below the threshold, all transverse 
perturbations are suppressed and TI is eliminated altogether. This is core idea of our Letter. 
The idea of using the threshold to eliminate TI applies to any type of nonlinearity, yet for the 
rest of this Letter we concentrate on the Kerr nonlinearity for two main reasons. First, wave 
propagation in Kerr media is described by the cubic nonlinear Schrodinger equation (NLSE) 
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which is one of the most general soliton equations [23]: it describes solitons (at least in some 
limit) in capillary waves, gravity waves, deep water waves, shallow water waves, plasma 
waves, and optics. Generally speaking, the NLSE describes envelope solitons in dispersive 
wave systems with weak symmetric anharmonicity. The second reason is that the effect of TI 
for Kerr solitons is very strong and we can demonstrate a convincing difference between 
having TI and eliminating TI by making the soliton incoherent along y. 

To formalize the problem, we use the radiation transfer method [17]. An incoherent beam 
can be considered as a series of coherent speckles that change, on average, every coherence 

time, icoh- We define B(x h y 1 x 2 ,y 2 , z) = <E (x 2 ,y 2 ,z,t)E(x h y 1 ,z,t)>, the spatial 
correlation function, where E(x,y,z,t) is the slowly varying amplitude. The brackets denote 
averaging over the response time of the medium, x, which is much larger than x CO h . From the 
paraxial wave equation we get [17] 

dB i ( d 2 B tfB , ik < x 

T~l ^T~k — + TIT" * = ~~ W x i ^ z)-ai(x 2 ,y 2 ,z)\B, ( 1 ) 

oz k oxop x oyopy n 

where z is the propagation direction, k the carrier wavenumber, n the bias refractive 
index, 8n the nonlinear contribution to the refractive index, x = {x 1 +x 2 )l2, y = (yi+y 2 )/2 the 
middle point coordinates, and p x =x r x 2 , Py=yi-y2 the difference coordinates. When p x =py =0, 
B(x,y,p x ,py,z) is the time- averaged light intensity, I(x,y,z). Let B s (x,y,p x ,Py) 

* 

= u(x + p x 1 2)u (x - p x I 2)4j (Py) be a stationary (z-independent) solution of Eq. (1). It 

represents a ID soliton stripe, which is self trapped and fully coherent in x, while being 
uniform and incoherent along y with an angular spectrum of A„(py). The waveform u(x) is 
determined by the nonlinearity, and can be taken to be real without loss of generality. In Kerr 
media, u(x) is proportional to sech(x). To study TI, we add a small perturbation Bj to B s 
where Bi«B s . The nonlinear index change in Kerr media is dn(I) = y I, where y is the 
nonlinear coefficient (n 2 ). Linearizing Eq. (1) yields 
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We seek solutions in the form B 1 (x,y,p x ,p y ,z)=exp(gz)exp(iay)L(x,p x )Af(py) + 
exp(g*z)exp(-iay) L*(x,-p x )A/(-p,), where a is the transverse wavenumber, g the TI growth 
rate (gain), and A/p,) the angular spectrum. The necessary condition B 1 (x,y,p x ,p y ,z) = 
B]*(x,y,-p x ,-py,z) is satisfied [17,19]. Substituting B } into Eq. (2) gives: 



gL(x, p x )A f (p ) - -{ v Hx ' A f (p ) + la 1 L(x, p x )} 
k dxdp x dp y 



ik 



[L(x + ^-,0)A f (0)e aPyl2 - L(x -^-,0)A f (0)e~' aPy ' 2 ]u(x + y)«(x - ^-)A n (p y ) 
+ [u(x + p x 1 2) 2 - u(x - p x 1 2f ]L(x, p x )A f (p y )A„ (0) 



(3) 



We are interested in determining the threshold condition, that is, to find the conditions 
under which the growth rate g(a) goes from a positive value to a negative value for all a. 
The procedure of determining the threshold applies to any form of spatial coherence (angular 
power spectrum), but for simplicity, we consider an initial Gaussian angular power spectrum, 

i.e. A n (p y ) - Qxp[-(p y k / 2) 2 ], where o defines the degree of coherence (the correlation 

distance). The higher o the more incoherent the soliton is. For a fully coherent soliton, if 
we were to calculate the growth rate g as a function of transverse wavenumber «, then g 
starts from (at a=0), increases and reach a maximum positive value (at the wavenumber 
with the largest growth rate), and then drops back to at a associated with the "cutoff 
wavelength" [2-4]. This means that for a coherent soliton, the growth rate is positive (and TI 
exists) for a band of wavenumbers a between zero and the cutoff wavenumber. For a soliton 
that is partially-coherent in y, we expect that if TI occurs, then g(a) will be positive for a 
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band of wavenumbers, just like in the coherent case. On the other hand, if TI does not occur, 
then g( a) will be negative for all a. In other words, if we start from G a small enough (a beam 
that is coherent enough), g(a) will be positive in a band of wavenumbers, just like the 
coherent case. But, as G a increases, this band becomes narrower until it completely disappears 
at some value G d t- If G is larger than this value, then transverse instability is completely 
eliminated. [Recall that modulation instability of incoherent light exhibits a similar trend: see 
Fig. 1 in [19], which shows g(a) of the incoherent MI as the threshold is approached]. We 
therefore expect that, at the threshold G a = G t, the two boundary points at which g(a)=0 
(one at a=0 and the other at the cutoff wavenumber), coincide. Thus, we seek the value of 
Go at which (I) g(a =0)=0 and (II) g'(a =0) = 0. We solve Eq. (3) by expansion while 
retaining up to 0(a). This eliminates the first term in the LHS because g = 0(a 2 ). For Kerr 
solitons, the growth rate of the transverse instabilities, g, is independent of x [2-4], that is, the 
growth rate of the transverse perturbations is the same across the ID soliton (even though 
the actual shape of the perturbations depends on x). We can therefore seek solutions of Eq. 3 
under the conditions (I) and (II) at the center of the soliton, i.e., at x =p x = 0, and safely 
assume that the threshold we find is the same everywhere on the soliton. It can be easily 
shown (by expanding into derivatives with respect to X; and x 2 ) that the second term of the 
LHS is zero for Xi = x 2 i.e. for p x = 0. Thus, from Eq. (3) we get: 



where I = u(0) is the peak intensity of the soliton. On physical grounds it is unlikely that a 
small perturbation will alter the coherence statistics of the soliton (especially here that 
propagation effects, given by g, are of order of a 2 and are ignored). Thus, we can safely 
assume that, A/py) = A„(py). Eq. (4) gives the threshold condition Gl T -2An m , where 
An = ]d is the maximum index of refraction change. One can actually calculate, using 

numerical methods similar to those of [10], the entire function g(a) and from it obtain the 
threshold for any angular distribution function [25]. 



k dp y 



Z(0,0) = °- yL(0, 0)ap y A f (0)exp[-(kG OT p y /2) 2 ] 



(4) 
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To verify the analytic predictions, we perform numerical simulations using the 
coherent density approach [12]. We launch a ID Kerr soliton with a Gaussian angular power 

spectrum, A(p y ) = exp[-(p y k 1 2) 2 ], for various values of . In this example, n =2.3, X 

in vacuum is 0.5 urn, FWHM = 9 urn, which yields a An = 0.0001056 and an analytic 
prediction of the threshold of 0T = 0.55°. Our numerical results are displayed in Fig. 1, 

where we show images of the intensity distribution of the soliton and cross sections of the 
intensity along y for x = 0. In Fig. 1(a) we show the input soliton at z = 0. In 1(b) we show a 
fully coherent soliton for (0 O = 0°) after 0.8 cm of propagation. As clearly depicted in the 
figure, TI breaks the soliton up into a train of 2D filaments. As we approach the threshold, 
the TI gain is getting smaller: as we set# to 0.5°, it takes a 4.5cm propagation to exhibit 
signs of TI [Fig. 1(c)]. To show that TI is completely eliminated when the nonlinearity is 
below threshold, we increase to 0.56° . As shown in 1(d), after 4.5 cm of propagation 
there are absolutely no signs of TI. 

Our experiments are conducted in photorefractive SBN:75 crystal in a setup similar to 
that of [12]. The beam is made spatially incoherent by passing it through a rotating diffuser. 
The rotating diffuser provides a new phase and amplitude distribution every x CO h ~ 1 |Lis, 
which is much shorter than the response time of the medium x ~ 1 s. Unlike all previous 
experiments with incoherent solitons, here we need to generate a beam which is very narrow 
and fully coherent in x, yet uniform and partially incoherent y. To do that, we use a cylindrical 
lens which focuses the beam only in the y direction onto the rotating diffuser. Then, by 
moving the focal point of this lens closer (farther away) from the diffuser, we increase 
(decrease) the coherence in y. The x coherence is not affected by the translation of this lens. 
After the diffuser, the beam is collimated (to ~ 2 cm) and passed through a narrow (along x) 
slit. The slit is made narrower than the speckle size in x, and it effectively creates a ID beam 
that is narrow and coherent in x and "infinitely" long (uniform) and incoherent in^. The slit is 
then imaged to the input face of the crystal. We get a reasonable estimate of the correlation 
distance by stopping the diffuser and measuring the average speckle size at the crystal input 
plane. Finally, we use an orthogonally polarized background beam that covers the crystal 
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uniformly as necessary for photorefractive screening solitons [8]. The input and output faces 
of the crystal are imaged onto a CCD camera. 

The photorefractive screening nonlinearity is in general saturable but resembles the 
Kerr nonlinearity when the peak intensity of the soliton normalized to the background 
intensity (or the dark irradiance) is much smaller than unity [8,24]. In our case, this ratio is 
~ 0. 1. At this normalized intensity, a soliton that is fully coherent in both x andy exhibits very 
strong transverse instability, as shown experimentally in [8]. We can therefore test the main 
prediction of our theory by gradually increasing the incoherence in y (decreasing the speckle 
size) until the soliton becomes transversely stable. All other parameters (applied field, 
intensity) are kept constant. Our results are shown in Fig. 2. The 12 um FWHM input beam 
[Fig. 2(a)] linearly-diffracts to a 60 um FWHM output beam after 6 mm of propagation in the 
crystal [Fig. 2(b)]. The nonlinearity is turned on with the application of 2.7 kV/cm and the 
beam self-traps forming a soliton in x. When the beam is fully coherent, then the soliton 
suffers from TI and breaks up into filaments [Fig. 2(c)]. When the beam is made incoherent in 
y, but with a large speckle size (-100 um), then the nonlinearity is still above the TI threshold 
and beam still suffers from TI [Fig. 2(d)]. Finally, by decreasing the speckle sizes to ~ 5 um, 
TI is eliminated and we get a stable (1+1)D soliton [Fig. 2(e)]. Thus, we have shown that a 
stable (1+1)D soliton can propagate in a 3D bulk medium provided that it is made 
partially-incoherent in the "uniform" transverse dimension. For the soliton to be stable, the 
degree of coherence in the "dimension of uniformity" must be such that the nonlinearity is 
below the threshold for transverse instability. 

In conclusion, we have derived theoretically the threshold for TI of (1+1)D solitons 
that are fully coherent in their direction of trapping yet are partially-incoherent in their 
direction of uniformity. We predicted that if the nonlinearity is below a well-defined 
threshold, then transverse instability of such ID solitons is completely eliminated. We proved 
our results analytically, numerically, and experimentally, and specifically showed that it is 
possible to generate stable ID Kerr-like solitons in a 3D bulk material. This is the only 
method we know of for propagating truly stable ID solitons in a bulk material. We note that 
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in principle our method applies to all types of saturable nonlinearities, and could be used to 
eliminate TI in them as well. Finally, we believe that this work opens up a whole range of 
possibilities of eliminating transverse instabilities in many soliton systems, for example, 
instabilities of ID dark solitons in bulk media, instabilities of ring beams (with and without 
topological charge) in self-focusing media, and more. 
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FIGURE CAPTIONS: 



Figure 1: Numerical simulations of a ID Kerr soliton with a Gaussian angular spectrum, for 
various degrees of coherence. Here, X= 0.5 urn, n =2.3, the beam FWHM is 9 urn, and An 
= 0.0001056 The analytically-predicted threshold is # =0.549°. (a) The input intensity, (b) 
The output beam after 0.8 cm of propagation for a fully coherent beam# = 0°, showing that 
the soliton is destroyed by TI. (c) The output beam close to the threshold (for = 0.5°), 
after 4.5 cm of propagation: as the threshold is approached, the TI gain is reduced and it 
takes a longer propagation distance before TI is evident, (d) The output beam for = 
0.56°, which is below the threshold, after 4.5 cm of propagation. TI is completely 
eliminated. 

Figure 2: Experimental results in a photorefractive SBN:75 crystal in the Kerr regime 
(intensity ratio ~ 0.1). (a) The input 12 urn FWHM beam, (b) The output beam after 6 mm 
of linear diffraction when nonlinearity is off. (c), (d), (e) show the output beams with an 
applied field of 2.7 kV/cm, for various degree of coherence along y and all other parameters 
kept constant, (c) A fully- coherent soliton breaks up into filaments because of TI. (d) The 
soliton is made incoherent along y but with large speckle sizes -100 um (small ) and still 
displays a strong TI. (e) The speckle sizes are reduced to ~ 10 um: TI is completely 
eliminated. 
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Fig. 1 

Charalambos Anastassiou et.al. 
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Fig. 2 

Charalambos Anastassiou et. al. 



